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Uniaxial elastomers are characterized by five elastic constants. If their elastic modulus C5 describing 
the energy of shear strains in planes containing the anisotropy axis vanishes, they are said to be 
soft. In spatial dimensions d less than or equal to 3, soft elastomers exhibit anomalous elasticity 
with certain length-scale dependent bending moduli that diverge and shear moduli that vanish at 
large length-scales. Using renormalized field theory at d = 3 and to first order in £ = 3 — d, we 
calculate critical exponents and other properties characterizing the anomalous elasticity of two soft 
systems: (i) nematic elastomers in which softness is a manifestation of a Goldstone mode induced 
by the spontaneous symmetry breaking associated with a transition from an isotropic state to a 
nematic state and (ii) a particular version of what we call a critically soft elastomer in which C5 = 
corresponds to a critical point terminating the stability regime of a uniaxial elastomer with C5 > 0. 

PACS numbers; 61.41. -|-e, 64.60.Fr, 64.60.Ak 



I. INTRODUCTION 

Liquid-crystalline elastomers Q, |^ |^ are elastic me- 
dia with the macroscopic symmetry properties of liq- 
uid crystals 4, 5]. They consist of weakly crosslinked 
polymeric networks with mesogenic units. The existence 
of the rubbery crosslinked network has apparently lit- 
tle impact on liquid crystalline phase behavior. In fact, 
the usual thermotropic liquid crystal phases, i.e., the 
nematic, cholesteric, smectic-A, and smectic-C phases 
have their elastomeric counterparts 0, @. However, 
because liquid-crystalline elastomers cannot flow, they 
have mechanical properties that differ significantly from 
standard liquid crystals. Usually, liquid-crystalline elas- 
tomers are prepared by crosslinking side-chain _,1Q] or 
main-chain |7| polymers. Alternative methods of syn- 
thesis include the polymerization of monomeric solutes 
in a liquid-crystalline solvent or the confinement of 
a conventional liquid crystal in a dilute flexible matrix 
such as aerosil 9j. 

The main subject of this paper are nematic liquid- 
crystalline elastomers, or briefly, nematic elastomers 
(NEs). For recent reviews on NEs see Refs. 0, 
\U\ . These materials have unique properties that make 
them candidates for device applications. Temperature 
change (12] or illumination |l3j | can alter the orientational 
order and cause the elastomer to extend or contract as 
much as 400% This qualifies nematic elastomers as 
contestants for use in artificial muscles ^1^, J^] . Another 
striking property of nematic elastomers is their soft elas- 
ticity [ij) ll^ I20I I2H characterized by vanishing shear 
stresses for a range of longitudinal strains applied per- 
pendicular to the uniaxial direction. 

The origin of soft elasticity in NEs is the spontaneous 
breaking of rotational symmetry of the isotropic state 
induced by the development of orientational order in 
the nematic state 0, Il8j | . This spontaneous symme- 



try breaking has the consequence that NEs are not like 
conventional uniaxial elastomers characterized by five in- 
dependent shear moduli. Rather, the elastic constant as- 
sociated with shear in planes containing the anisotropy 
axis (customarily called C5) vanishes in NEs. 

Though we are primarily interested in NEs, we also 
study a particular version of another class of soft uni- 
axial elastomers in which the elastic constant C5 simply 
vanishes for energetic or entropic reasons. In this case, 
C5 = is a critical point marking the boundary between 
the high-symmetry uniaxial phase with C5 > and a low- 
symmetry sheared phase with C5 < 0. In other words, 
C5 acts like the temperature variable in a standard ther- 
mal phase transition. A complete model of this critical 
point requires the introduction of third and fourth or- 
der terms in the nonlinear strains to stabilize the system 
when C5 < 0. Since this complete model is character- 
ized by a large number of parameters, it is quite com- 
plex. Rather than analyze this full model, we consider 
simpler model systems, which we call critically soft elas- 
tomers or CSEs, defined by simply setting C5 = in the 
standard elastic energy of a uniaxial medium containing 
only quadratic terms in nonlinear strains. Remarkably, 
CSEs exhibit well-defined anomalous elasticity much like 
that of the more physical NEs even though they lack the 
nonlinear terms needed to stabilize their low-symmetry 
phase. CSEs are simpler in many ways than NEs, and 
the analysis of their anomalous behavior provides a useful 
and instructive tutorial prelude to the analysis of NEs. 

On the level of mean-field theory, the elastic energies 
of NEs and CSEs coincide. Therefore, their respective 
elastic properties are equivalent above their mutual up- 
per critical dimension 3. For dimensions d < 3, however, 
fluctuations become important. These fluctuations lead 
to Grinstein-Pelcovits |22l | type renormalizations culmi- 
nating in anomalous elasticity, i.e., in a length-scale de- 
pendence of certain elastic constants, with different uni- 
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versality classes for NEs and CSEs. 

In the present paper we explore the anomalous elastic- 
ity of NEs and CSEs by carrying out a renormalization 
group (RG) analysis. Using the methods of renormalized 
field theory we examine the scaling behavior of the 
elastic constants in d = 3 as well as in d = 3 — e dimen- 
sions. A brief account of our work on NEs appears in 
Ref. 24] . Our work on CSEs has not been reported hith- 
erto. We will treat only systems in which random stresses 
are not important. Random stresses lead to a different 
universality class with anomalous behavior [r^l23 | below 
five rather than three dimensions. 

The plan of presentation is as follows: In Sec. |n] we 
give a brief summary of our main results. In Sec. IIIII 
we discuss CSEs. In Sec Him we briefly review some 
elements of the Lagrangian theory of elasticity and then 
set up a Landau- Ginzburg- Wilson elastic energy func- 
tional (Hamiltonian) for CSEs. As a prelude to the sub- 
sequent RG analysis we analyze the symmetry contents 
of this Hamiltonian. Section IIII Bl contains the core of 
our renormalized field theory of CSEs. We explain our 
diagrammatic perturbation calculation and its renormal- 
ization. By solving the appropriate RG equation, we as- 
certain the scaling behavior of displacement correlation 
functions and ultimately that of the relevant elastic mod- 
uli. We conclude Sec. lIIII bv making contact with conven- 
tional uniaxial elastomers by incorporating a small but 
nonvanishing C5, which leads to semi-soft elasticity. Sec- 
tion llVl deals with NEs and has an outline similar to that 
of Sec. IIIII In Sec. IIV Al we derive a Landau-Ginzburg- 
Wilson minimal model for NEs in form of a field theo- 
retic Hamiltonian. Our renormalized field theory for this 
model is presented in Sec. IIVBI The main part of this 
paper concludes with Sec. where we give concluding 
remarks. There are five Appendices. In Appendices 1X1 
lEl and 10 we derive Ward identities for CSEs and NEs. 
An alternative RG approach to NEs is sketched in Ap- 
pendix^ Appendix^contains details on the calculation 
of Feynman diagrams. 



and Udd, Uab, and Uad are c omp onents of the Lagrangian 
nonlinear strain tensor |26l l27l | . The elastic constant Ci 
describes longitudinal shear along the anisotropy. C2 
couples strains along the anisotropy axis to shears in the 
plane perpendicular to the anisotropy axis. The elastic 
constants C3 and C4 are associated with shear purely in 
the plane perpendicular to the distinguished direction. 
K is a bending modulus. The pure CSE system with 
soft elasticity is characterized by C5 = 0. The coupling 
C5 plays a role similar to that of the temperature in a 
thermal phase transition. When C5 > 0, the system 
displays conventional uniaxial elasticity. When C5 < 0, 
the system is unstable with respect to the formation of 
a lower-symmetry (sheared) elastic state. C5 = marks 
the transition between the two phases. At this point, the 
system shows critical behavior analogous to critical be- 
havior at a thermal phase transition. Our simple CSE 
model does not include higher order terms in the strain 
to stabilize the sheared phase. We expect, however, that 
the full model with these terms included will have the 
same structure as our CSE model, at least for C5 > 0. 

The inverse displacement correlation functions (vertex 
functions) obey scaling forms that can be expressed for 
d < 3 as 
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II. SUMMARY OF RESULTS 

For the convenience of the reader we now summarize 
our main results before we get into details of our work. 

A. Critically soft elastomers 

The elastic constants of CSEs are defined via the model 
elastic energy 

H=^J df'^x^ j dxd {Ci + K {Vludf (2.1) 

-t- 2C2 UddUaa + C3 + 2C4 UabUab + C'5 UadUad^ , 

where Ud and MQ,a=l,...,d— 1, are, respectively, the di- 
rections parallel and perpendicular to the nematic order 



where qd and are, respectively, wavenumbers (or mo- 
menta) parallel and perpendicular to the anisotropy axis. 



7JK = 4£/7 , 77c = e/7 , TJ2 = e/7 , (2.3a) 
1^5 = 1/2-6/28, (4-77k)/2, (2.3b) 

where £ = 3 — d, are scaling exponents. In principle, 
the exponents (p and 772 could be independent of the 
other exponents. For the CSE model, however, they are 
not at least to first order in e. The scaling functions 
^dd{ci±,qd,e), $arf(q_L,gd, e), and l>ab(q_L, gd, e) are re- 
spectively proportional to q^ -I- g^ -I- e q^ , q^gd and q^ 
in the long-wavelength limit in mean-field theory. $afc 
also has a term proportional to g^, but its coefficient 
is irrelevant and we will not be concerned with it here. 
The boundary between scaling and Gaussian behavior is 
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marked by the non-linear length scales 
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where T is the temperature measured in units so that the 
Boltzmann constant is equal to 1. 

The above scaling forms imply the following scaling for 
the elastic constants: 



C-2 
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^ C4 (2.5) 

{Ldqd)'''^'* if 6|q±|irf/ii » 1 , qj. = 
[L^i-^^c if (q^,g,) = (O,0) 



(L^lq^l)-"- if C5|q±|»l, gd = 
(idgd)-"^/* if Cslq^l^d/ii » 1 , = 
(i±4"')~"" if (qi,g<i) = (O,0) 

(2.6) 



where ^5 is a correlation length given by 



(2.7) 



The elastic constant Ci is not renormalized and it is not 
singular in either the wavenumbers or ^5. When C5 is 
nonzero, T^d ~ C2^q\ + Kq\ at small q, where 75 — 
(2 — riK)v5 and K is given by the last expression in Eq. 
(EH)- 

At exactly three dimensions the above power law sin- 
gularities become logarithmic singularities 

C2 ~ C3 ~ C4 ~ I Indq^l^)!-!/^ , (2.8a) 
^ |ln(|q^|^)|V^ (2.8b) 

where /i is a wavenumber scale. This logarithmic 
anomaly can be observed for ^±|qj^| ^ 1, where 
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exp 



(2.9) 



In the critical regime at small C5 , we find the following 
universal Poisson ratios: 

Ci/{CiC4) = and C3/C4 = -l/2. (2.10) 



B. Nematic elastomers 

At the transition from the isotropic to the ne- 
matic phase, liquid-crystalline elastomers undergo an 
anisotropic stretch relative to their isotropic reference 
state of a factor Aqu along the nematic axis and Aq^ 
perpendicular to it. In the absence of explicit uniaxial 
terms such as 
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(2.11) 



the elastic energy of the nematic phase that forms spon- 
taneously from an isotropic phase is soft, i.e., the elastic 
constant C5 for shears in the uniaxial plane vanishes. Af- 
ter rescaling lengths measured relative to the anisotropic 
nematic reference state and displacements according to 

Xa Xa, Xd —> \/r - IXd, Ua -> Ua and Ud \/r -1 Ud, 

where r = A^h/Aqj^, the relevant parts of the elastic en- 
ergy can be written when uniaxial terms are present and 
small as 

H = i Jd^^x^ j dxd {Civil + K {Vlud)\2.l2) 

+ 2 C2 VddVaa + C3 + 2 Ci vl, + C5 via} > 

with nonstandard strains Vab — ^{daUb + d^Ua — 

daUddbUd), Vdd = ddUd + \daUddbUd and Vad = \daUd 
and where C5 goes linearly to zero with h or the magni- 
tude of other uniaxial terms. 

The scaling of the inverse displacement correlation 
functions of NEs is similar to but not identical to that of 
the CSE model: 
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where, to first order in e 



TjK = 38e/59 , T]c = 4e/59 , (2.14a) 
(j) = 2- 21e/59 , iyh = 1/2 + 9e/108 , (2.14b) 



and 



Li_ 
Ld 
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(2.15a) 
(2.15b) 



Note that four independent scaling exponents, ijk, rjc, 
(f), and Vfi are required to describe NEs with a small uni- 
axial energy. In the above, lengths, displacements and 
q- vectors are measured in rescaled units. 

The above scaling forms predict that Ci, C2, and C3 
are unrenormalized and that 

(iilq^l)"^ if 6.|q±l»l, qd = 
Ci ~ i {Ldqd)''^'^ if ih\ciALd/Li_ » 1 , q^ = 
{L^i-^yc if (q^,5^)^(o,0) 

(2.16) 
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as well as 

(L^lq^l)-"- if a|q±l»l, qd = 
K ^ < {Laqd)-''"^'^ if a|q±IW^± » 1 , = 
(i^C')""" if (q^,gd) = (O,0) 

(2.17) 



with £^h given by 



a = {L\h/K)- 



(2.18) 



At small but nonzero V^^ ^ h?'^q\ + Kq\ at small 
q, where 7/1 = [2 — rjh)i'h and X is given by the last 
expression in Eq. (|2.17|) . 

At exactly three dimensions the above power laws be- 
come 



C4 ^ llndq^l/^)!-^/^^ 
^ |ln(|q^|/M)r/59. 



(2.19a) 
(2.19b) 



The length scale that marks the crossover from harmonic 
to logarithmic behavior is 



= ^ exp 



(2.20) 



Provided that C5 is small, the critical regime entails 
four independent Poisson ratios: 

C2/C1 = 1 , C3/C1 = 1 , C4/C1 = , (2.21a) 
(2C2 - C3 - Ci)/Ci = 1/2 . (2.21b) 

III. CRITICALLY SOFT ELASTOMERS 

A. The model 

We start by setting up a field theoretic minimal model 
for CSEs that is suitable for our subsequent RG anal- 
ysis. We find it convenient to use the Lagrangian for- 
mulation of elasticity "2^, "271. In this formulation the 
mass points of the equilibrium undistorted medium are 
labeled by their position vectors x in d-dimensional (ref- 
erence) space. When the medium is distorted, a mass 
point originally at x is mapped to a new point R(x) in 
d-dimensional (target) space. Since R(x) = x when there 
is no distortion, it is customary to introduce the phonon 
variable u(x) = R(x) — x that measures the deviation of 
R(x) from x. 

Suppose for a moment the medium is distorted solely 
by stretching. The energy of the distorted state relative 
to the reference state depends on the relative amount of 
stretching dR? — dx^ = 2uijdxidxj , where 



diUkdjUk} 



(3.1) 



with i,j,k — 1, . . . ,d [2^ are the components of the fa- 
miliar nonlinear Lagrangian strain tensor u. Note that 
u is invariant under arbitrary rotations in target space. 



This feature makes the Lagrangian strain tensor an ade- 
quate variable for formulating elastic energies because all 
elastic media are rotationally invariant in target space. 
This invariance is easy to understand: different physi- 
cal orientations of the same sample have the same en- 
ergy [23. As customary, the reference state, relative to 
which u is defined, is taken to be in mechanical equi- 
librium, and hence no terms linear in appear in the 
stretching energy. To lowest order, the stretching energy 
is then of the form 



^^st = ^ / d'^xKijkiUijUki 



(3.2) 



where Kijki is an elastic constant tensor. For media 
isotropic in the reference space, for example, there are 
only two independent elastic constants in isTyfej that are 
known as the Lame coefficients A and /i. Media with uni- 
axial symmetry in the reference space are characterized in 
general by five independent elastic constants. Assuming 
that the anisotropy axis is in the = (0, . . . , 1) direction 
we may write the stretching energy as 



d'^^X±_ J dXd{Ciu\^ + 2C2UddUaa 



(3.3) 



with d^ — d ~ \ and a,b = 1, . . . , c?j_. 

Now suppose that the elastic constant C5 vanishes. 
Rewriting Hgt in Fourier space one sees easily that the 
stretching energy cost is zero for phonon displacements 
u(q) perpendicular to with momentum q parallel to 
ed and for u(q) parallel to with q perpendicular to 
Bd- In other words: CSEs are soft elastic materials. 

For many elastic systems it is justified to neglect ener- 
getic contributions, such as bending, that are associated 
with higher derivatives of the displacements. That is, be- 
cause bending is unimportant compared to stretching at 
small momenta. Due to the soft elasticity, however, the 
stretching energy of CSEs can vanish, and hence, bending 
is important. 

For the moment, we set aside the uniaxial term propor- 
tional to C5 and concentrate on the pure soft case. The 
effects of a small but non- vanishing C5 will be included 
later on. Taking into account stretching and bending, 
the CSE model is defined by the Hamiltonian 

n^^Jd^^x^ Jdxd[Ci ul^ + K{vludy 

+ 2C2 UddUaa + C3 + 2C4 UabUabj , (3.4) 

where AT is a bending modulus. All other bending terms 
allowed by symmetry turn out to be irrelevant in the 
sense of the renormalization group. Also, not all parts 
of the strains are relevant. Discarding any parts of the 
strains that lead to irrelevant contributions to the Hamil- 
tonian, as discussed further below, leaves us with 
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Uab 



{daUb + dbUa + daUddbUd} 



(3.5a) 
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and 

Udd = ddUd ■ (3.5b) 

In principle we could use H as it stands in Eq. (|3.4|l 
for our RG analysis. We find it convenient, however, 
to reduce the number of constants featured in H at the 
onset. To this end, we rescale Ua {K /C/CjUa-, Ud 
\J KjC^Ud, and Xd — > {Ci/ K'^)xd- Then the Hamiltonian 
takes on the form 

'^ = \j d"^^^^ J dxd[ujula+ {Vludf 

+ 2gUddUaa + fula+^UabUabj , (3.6) 

where 

Lu = C\Kycl, g^C2{K/Cif/\ f^Cs/CU. 

(3.7) 

At this stage we would like to point out that we explic- 
itly keep the temperature T in the Boltzmann weight 
exp(— 7^/T) [30| governing our field theoretic calcula- 
tions. In what follows we carry out a perturbation expan- 
sion in the temperature, i.e., T serves as our expansion 
parameter. As a consequence, not only the constants and 
fields featured in H but also the temperature will require 
renormalization. 

An effective Hamiltonian for Ud alone can be obtained 
by integrating out the transverse variable Ua from the 
full CSE Hamiltonian of Eq. (O. When C4 oo, this 
process leads to the Hamiltonian for a smectic-A liquid 
crystal whose anomalous elasticity was analyzed by Grin- 
stein and Pelcovits 22] . Our rescaling of variables to ob- 
tain Eq. 1)3. 6|l with the coefficients of both (V^u^)^ and 
2uabUab sct to Unity is not ideally suited to taking the 
C4 = 00 limit. Our primary interest is the anomalous 
elasticity unique to soft uniaxial systems for which the 
parametrization of Eq. H3.6|) is appropriate. We will not 
give further consideration to the Grinstein-Pelcovits limit 
of our model. 

As a further step towards our RG analysis, we now 
discuss the scaling symmetries of our model. First, under 
a global rescaling of the coordinates Xa ^~^Xa and 
Xd ~* ^■^'^Xd, we find a scaling invariant theory provided 
that Ua — > /i^Ua and T — > fj.^^'^T (/i invariance). Viewing 
fj, as an inverse length scale, this means that the field Ua 
has a naive dimension 1 and that the naive dimension of 
T is e = 3— d. The field Ud and the remaining parameters 
in Ti. have a vanishing naive dimension. Above d = 3 
dimensions the naive dimension of T is negative and T 
is irrelevant whereas it is relevant below d = 3. Hence, 
we identify dc = 3 as the upper critical dimension of the 
CSE model. 

At this point we take a short detour and catch up 
on justifying the truncation of the strains as stated in 
Eqs. 1)3. 5|l . Applying the /x-rescaling to the original full 
strains leads to 

Uab {daUb + dbUa + daUdObUd + daUcObUc] 

(3.8a) 



and 

Udd S ddUd + — {ddUdf + —ddUcddUc > .(3.8b) 

The terms carrying extra powers of ^ do not contribute to 
the leading behavior in the limit /i — + 0. Hence, they can 
be neglected in studying the long length scale behavior 
at leading order, i.e., they are irrelevant in the sense of 
the RG. 

Second, due to the anisotropy of the model, we may 
rescale the longitudinal coordinate alone: Xd — > Pxd- 
Scale invariance is retained if ijj — > 0^uj, g ^ fig, and 
T f3T (/3 invariance). Note that the composed cou- 
plings 

cr = g'^/uj, p = f, and t = ^"^T/^^; (3.9) 

are invariant under the longitudinal rescaling. The fac- 
tor 11^^ is included in the definition of t to render it, 
like tJ and p, dimensionless. As we go along, we will 
see that tr, p, and t emerge quite naturally in pertur- 
bation theory. Third, Ti. is invariant under the rescal- 
ing Ud ^ Ud + fd{^±), where fd is an arbitrary func- 
tion of the transversal coordinate x^. Fourth, rescaling 
Ua Ua + fa{xd) + MabXb leavcs Ti. invariant if the ma- 
trix constituted by the Mab is antisymmetric and fa is 
a function of the longitudinal coordinate only. Finally, 
Ti. is invariant under the transformation Ua Ua + OdUd 
and Ud Ud — OaXa provided that the 0's are small. Note 
that this transformation mixes the longitudinal and the 
transversal fields (mixing invariance). It can be viewed 
as a remnant of the rotational invariance of the original 
theory in target space. This mixing transformation will 
be valuable for us because it leads to Ward identities that 
reduce the number of vertex functions to be calculated 
in perturbation theory. These Ward identities will be 
derived in Appendix IXI 



B. Renormalization group analysis 

In this section we determine the scaling behavior of the 
correlation function of the fields Mq(x) and Ud(x) by us- 
ing perturbation theory augmented by renormalization 
group methods. As usual, we analyze vertex functions 
that require renormalization due to the presence of ul- 
traviolet (UV) divergences in Feynman diagrams. Our 
main tools in this section will be dimensional regular- 
ization and minimal subtraction. To avoid infrared (IR) 
singularities in the Feynman diagrams, we supplement 
our Hamiltonian with a mass term, 

n^H+^ J d'^^xx_ j dxd ul . (3.10) 

At the appropriate stage of the calculations we then sent 
T to zero to recover the original situation. 
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1. Diagrammatic expansion 

In order to set up a diagrammatic perturbation ex- 
pansion we have to determine its constituting elements. 
First, we have the Gaussian propagator G that has the 
form of a (i X d matrix. The elements of its inverse F are 
readily collected from the Hamiltonian: 



F 



dd 



^Qd 



Tad = T gqaQd , 

Tafc = T-^[if + l)qaqb + Sabqii] 



(3.11a) 
(3.11b) 
(3.11c) 



Inverting F we find that the Gaussian propagator has the 
elements 



Gdd 
Gad 
Gab 



T- 



Bt 



Bt 



= T 



B 



Ml 



qaQd 



Aql + Bci\ q2 



6ab Dt + C q^ + D<i\ 

QaQb 



Bt + Aq^ 



Sql q1 



(3.12a) 
(3.12b) 
, (3.12c) 



where we have used the shorthand notations A = uj{f + 
2)~g^, B = f + 2,C ^ uj{f + 1) -g\ and = / + 
1. Second, our diagrammatic expansion features the 4 
vertices 



X'^X 



■ 9 (1) (2) (3) 
^^Id % % 

■ f (1) (2) (3) 

'Qb % 

■ 1 (2) (1) (3) 



T 



la lb : 



/ + 2 g(l)j(2) (3) (4) 
la Ha % lb 



(3.13a) 
(3.13b) 
(3.13c) 
(3.13d) 



It is understood that the sum of the momenta has to 
vanish at each vertex. 

Next, we need to determine which of the vertex func- 
tions F^^^'^^ with M external Ua-legs and N external 
Wrf-legs are superficially UV divergent. Analyzing their 
topology, we find that the superficial degree of divergence 
S of our diagrams is given at the upper critical dimen- 
sion by (5 = 4 - M - 2£)|| - D^, where Dy {Dj_) is the 
number of longitudinal (transversal ) derivatives on the 
external legs. Thus, the only vertex functions containing 



superficially divergent diagrams are T^a'^\ 



.(1,1) ^(2,0) 



ad 



ab 



p(0,3) p(l,2) 
ddd ' ^ . 



and F^"'^^ 
add ' '^'"^ ^ dddd- 



.(0,2) 
dd 



All these vertex 



functions have to be taken into account in the renormal- 
ization procedure. By virtue of the mixing invariance, 
however, there exist several relations between the vertex 
functions in form of Ward identities. These are derived 
and stated in Appendix IXI Due to these Ward identities 
it is sufficient for our purposes to actually calculate the 



two-point functions F 
equations of state F^"'^ 



(0,2) ^(1,1) 



dd 



ad 



and Fi^^°) 



and Tfi^'^K Once the 



are satisfied 

and the 2-point functions are renormalized, the Ward 



FIG. 1: Feynman diagrams contributing to Tdd- The dashed 
lines symbolize the element Gdd of the Gaussian propagator. 
Lines half dashed and half solid with an index, say a, stand 
for Gad- Solid lines with 2 indices, say a and b, visualize Gab- 
The ticks indicate derivatives with respect to the reference 
space coordinate with an index specifying the component. 



d 



I d 



FIG. 2: Feynman diagrams contributing to Vad- The meaning 
of the symbols is the same as in Fig. 



identities guarantee that the remaining vertex functions 
are cured of their UV divergences. 

We calculate the two-point vertex functions to one-loop 
order using dimensional regularization. The Feynman 
diagrams entering this calculation are listed in Figs.^to 
13 Details on computing the diagrams can be found in 
Appendix^ Our results for the 2-point functions read 
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FIG. 3: Feynman diagrams contributing to Tab- The meaning 
of the symbols is the same as in Fig. Q 



(3.14a) 



-e/4 



167re 

3(1 



i2 + f)L0 

-/)ql 



.(1,1) 



T 



V2- 



9 QaQd 

7(1 + /) 5' 



(3.14b) 



-e/4 



167re^-g^ + {2 + f)ij 



.(2,0) 



V2- 



T-'[{f + l)qaqb + 6ab qi 

7[2(l + /)%agb+5a6qi] _ 



(3.14c) 



327:e ^-g^ + {2 + f)uj 



e/4 



As already indicated by the vertex functions' superficial 
degree of divergence, higher order terms in the momen- 
tum expansion are convergent and hence can be neglected 
for our purposes. Likewise, contributions proportional 
to the mass r are convergent. Thus, r does not require 
renormalization and consequently its scaling dimension is 
identical to its naive dimension four. At this stage of the 
calculation r has fulfilled its purpose, viz. it prevented 
IR singularities from producing spurious e poles. In con- 
junction with the e expansion we now can safely send r 
to zero. 



inal structure: 



n 

T 



+ 2Zgg UddUaa + Zff ul^ + 2 ZTZUabUab^ , (3.16) 



The simplest way of determining the renormalization Z- 
factors is minimal subtraction. In this procedure the Z- 
factors are chosen so that they solely cancel the e poles 
and otherwise leave the vertex functions unchanged. Ex- 
pressed in terms of the effective couplings introduced in 
Sec, nil A1 our Z-factors are of the structure 



Z,„{t,a,p) = 1+51 



xM(t,a,p) 



(3.17) 



m— 1 



The A"'^™^(t, (T, p) are expansions in the effective temper- 
ature t beginning with the power t™. It is a fundamen- 
tal fact of renormalization theory, cf . Ref. 23] , that this 
procedure is suitable to eliminate all the UV-divergences 
(not only the superficial ones) from any vertex function 
order by order in perturbation theory. To one-loop order 
we find that our Z factors are given by 



Z = l + t 



ZT = l-t 



14+ 13p 



32 TT £ ^/2Tp V2 + p - cr 

3(1 + P) 
\^Te^/2Tp^/2 + p- a 



Z^ = l + t 



16 7r£ ^/2 + p — cr ' 



z, = i+t^/I±pMA, 

" SlTey/2 + p-(J 



Zf = l + t 



x/2Tp(l + 4p + 2p^) 
32TTe^2 + p-a 



(3.18a) 
(3.18b) 
(3.18c) 
(3.18d) 
(3.18e) 



2. Renormalization 



The UV divergences have their manifestation in the e 
poles appearing in Eqs. (|3.14() . We eliminate these poles 
by employing the renormalization scheme 

Ud ^ Ud ^ Z^/^Ud, (3.15a) 

Ua ^ Ua = ZUa , (3.15b) 

T^f = ZZj^^T, (3.15c) 

uj ^ uj — Z^^ Zi^io , (3.15d) 

g^ g^ Z-'^^Zj^'Zgg, (3.15e) 

f-^f = Z-'Z^'Zff, (3.15f) 

where the ° indicates unrenormalized quantities. Our 
scheme is chosen so that the Hamiltonian retains its orig- 



3. Scaling 1: RG equation and its solution 

Next, we infer the scaling behavior of vertex function 
from a RG equation. This RG equation is a manifesta- 
tion of the fact that the unrenormalized theory has to be 
independent of the arbitrary length scale introduced 
by renormalization. By virtue of this independence, the 
unrenormalized vertex functions satisfy the identity 

^^t(M,N) ({q^, ^4 ; ti, f , ff, /) = . (3.19) 
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The identity (|3.19() translates via the Wilson func- 
tions 



^ 1 7 

In Z,. 



dlriL 



It - lu. 



1 



1 



Pa = fJ. 



da 
dp 



7 + + -^lu 



P {l + lT-lf) 



(3.20a) 
(3.20b) 
, (3.20c) 
(3.20d) 
(3.20e) 



into the Gell-Mann-Low RG equation 



D„ 



M + - ) 7 



(3.21) 



Here we have used the shorthand notation 



d 



d 



d 



d 



d 



The Wilson-7-functions are easily gathered from the 
renormalization factors stated in Eqs. H3.18|l upon re- 



expressing p-^ as /3) 



'* at ■ 



Because the Wilson functions 



must be finite, one then immediately gets 

7...(i,fT,p) - -tatX«(t,(7,p), (3.23) 

where X'^'^^ is defined in Eq. H3.17|l . Since we will need 
them to determine the fixed points of the RG flow, we 
state the Wilson-/3-functions explicitly: 

20 + 19p -2a- pa 



Pt = -te + r 



32 7ry2+7V2 + p- CT 



Pa = t 
Pp = t 



V2Tpcr(3 + 4p-2g) 
327r V2 + p - CT 
2 + 7p + lp^ + 2p3 



, (3.24a) 
(3.24b) 
(3.24c) 



m-K ^/2 + p ^/2 + p - a 

To solve the RG equation we employ the method of 
characteristics. We introduce a flow parameter ^ and 
look for functions p{t), Z{£), u){t), t{i), a{e), and p{£) 
determined by the characteristic equations 



dl 



p, p{l)^p, 



(3.25a) 



e^^\nZ{e)^jm,a{£),p{l)) , Z(l) = l, (3.25b) 
i^^luQie)^ a{l),p{e)) , c^(l) = , (3.25c) 



(3.25d) 



ig-^m = Pa iaii),m,m) , '^(l) = ^ , (3.25e) 

i§-gm = pp m,m,m) , m = p ■ (3.25f) 



These characteristics describe how the parameters trans- 
form if we change the momentum scale p according to 
p — > p{£) = £p. Being interested in the IR behavior of 
the theory, we focus on the limit £ ^ 0. In this IR limit 
we find that the set of coupling constants {t{£), a{£),j5{£)) 
fiows to a stable fixed point 



/32 7re 



0,- 



(3.26) 



Pa{t*,a*,p*) 



[t ,a ,p ) 

satisfying Pt{t*,a*,p*) 
Pp{t*,a*,p*) = 0. Recalling that a = g'^/ui and 
p — f, we learn that the stable fixed point implies two 
universal ratios of the elastic moduli (Poisson ratios): 

Cl/{CiC4)^0 and €3/04 = -1/2. (3.27) 

Note that the longitudinal and the transversal directions 
are effectively decoupled at the IR stable fixed point. In 
addition to the stable fixed point there are 4 unstable 
fixed points, viz. the zero temperature fixed point t* = 
as wefl as (327re, 0, -1), (64^2/3 7r£/13, 1/2, -1/2), and 
(32^273716, -1/2,-1). 

With help of the characteristics the RGE is readily 
solved, at least formally: 

r^^-^'^) i{ci^,q,};^A'y,P,p)^Zi£)-^^'+^/'^ 
X r(*^'^) ({q±, <z4 ; ioi£),m,^{^), m,pi) -(3.28) 

As it stands, Eq. H3.28|l does not account for the naive 
dimensions of its ingredients. Recalling the p invariance 
of the Hamiltonian H it is straightforward to check that 

r(*^-^) {{q^,q,};cu,t,a,p,p) 

-(d+l)+dM+{d+l)N 



= M 



— } ■,Lj,t,a,p,l 



q 

p p 

Moreover, the P invariance tells us that 
r(^'^'^) {{q^,q,};LU,t,a,p,p) 



(3.29) 



(3.30) 



- --(^'^+^-i)/2r(*^^^) [{q^,u;'/\,}-l,t,a,p,p) . 

where the last line reflects our freedom to chose P = 
uj-^l"^. Combining Eqs. ijX^ . (jX^ we find that the 
scaling behavior of the vertex functions is described by 

r^*"'"^) ({q^,g4;^,i,^,p,A^) 

Qd 



[\p£' 



(3.31) 

■,u;{£),t{£),m,P{i),pe^ 



{p£Y 

Due to Eq. (|3.30() we may also write 
r(*'''^)({q±,gd};w,t,a,p,A*) (3.32) 



{p£f 



:Xm,a{£)rp{t),pi 
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We have to mention that Eqs. (|3.31() and (|3.32() do not 

correctly describe r|j^'°-' at qj^ = because we omitted a 
dangerous irrelevant bending of the type Kd q^- Without 
such a term, Eqs. H3.31|l and H3.32|l unphysically suggest 
that the leading scaling behavior of r^j^'"' is independent 
of qd for q^ =0. In other words, Kd is dangerous irrel- 
evant as far as the leading behavior of v'^^^f'^ at q^ = 
is concerned. Because Kd is irrelevant its omission has 
no impact on the leading behavior of the relevant elas- 
tic constants and hence does not affect our main results. 
Of course one could and it would be interesting to inves- 
tigate the scaling behavior of Kd- One has to keep in 
mind, however, that irrelevant terms tend to mix under 
renormalization with a whole bunch of other irrelevant 
terms, making a proper RG analysis a tedious endevour. 
This will be left to future work. 



4- Scaling 2: Physical quantities 

The variables we considered so far in our RG analysis 
had the benefit of being convenient. The flip side of this 
convenience is, however, that the featured quantities have 
no direct physical meaning. Now we recast our results so 
that their physical content becomes pronounced. 

To have a clear distinction between physical variables 
and the scaled variables we used in our calculations we 
mark the latter in the remainder of this section with a 
hat, i.e., we denote 



Id 



K^ 



qd, —U. 

K 



Ud 



—Ud- (3.33) 



It is not difficult to see that the relation between the 
physical vertex functions and the vertex functions in the 
scaled variables is given by 

T^'''''H{ci±,qd};u;,t,a,p,f,) (3.34) 

^ ^M+3A./2-2c.l-A'/2f(AW ^{q^, q,} , ^ a, p, f,) . 



Blending Eq. H3.34|l and our findings of Sec. IIII B3l we 
now obtain 

r(*^^^) ({q^, g4 ; CO, t, a, p, M) = i c;^M^^-'^/2 

^ ^{2M+N~2)/2-^^(M~2)/2j^2(d^l)~dM~{d+l)N 
^ ^-{d+l)+dM+{d+l)N ((:^(£)/t^)-(^^+^-l)/2 ^(•^)-(M+W/2) 
(M,N)( \ L^Ci±_ ^Ji^(l)I^Ldqd 



X (f*^ 

l,t(^),a(^),p(£),lj , (3.35) 
where we introduced the susceptibilities 




to make explicit that the vertex functions are propor- 
tional to the inverse temperature. Moreover, we intro- 
duced the length scales 

(3.37a) 
(3.37b) 

At this stage, L_l and Ld are still arbitrary. Further be- 
low we will fix these length scales so that they acquire a 
physical meaning, viz. the borderline between harmonic 
and scaling behavior. Note that the result H3.35|l is gen- 
eral in the sense that it holds in the harmonic as well as 
in the scaling limit. When t(X) 0, the system behaves 
approximately like a harmonic system. When I is small, 
behavior is determined by the fixed point with t* ^ t. 

Now to our main goal, viz. the behavior of the elas- 
tic constants. The sought after behavior can be inferred 
without much effort from result (|3.35(l . As an example, 
we consider the case M = 0, = 2 in some detail. In 
addition to the information contained in Eq. H3.35() we 
need some knowledge on the concrete form of the scaling 



function <1> 



(0.2) 
dd 



Since t* is of order e it is reasonable to 



assume that '^'^df' '-'^^ approximated by its Gaussian 
form [cf. Eq. H3.11|l ] even in the scaling limit. Hence, we 
write 



(0,2 
dd 



^({^, ^^^-'^ -^};M-(^),^(^),p(^),: 

■ ^_ \2 ^ 

^jL0{l)/ujLdqd\ /i_Lq_L\ 



(3.38) 

Merging Eqs. (|3.38|) and H3.35|l we obtain the physical 
vertex function 

ri°/^ {Ci^.qd) = {C,ii)qj+Kii) |q^|4} (3.39) 



with 



K{£) = Kr'Z{£)-^ [t{i)/t] [<^(^)H 



-1/2 



(3.40a) 



1/2 



(3.40b) 



In the case of the bending modulus we cannot finalize 
our conclusions without solving the characteristics, i.e. 
without knowing Z{£) and so on. For Ci, however, we 
observe without further information the following: 



4 CM 



7- 



t 



-7a. 



(3.41) 



Since 7^ is proportional to a it vanishes at the stable fixed 
point. Hence, Ci is independent of ^. In other words, Ci 
is normal. Of course, we cannot tell from our analysis 
if this stays true beyond 1-loop order. Nevertheless, this 
may well be the case. 
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The behavior of the remaining elastic constants can be 
extracted by similar means from the other 2-leg vertex 
functions. We find 

ri'/^ (q±, -Zrf) = T-^C2{i) qaqd , (3.42a) 

+ C4(^)(ga% + '5„fc|qj_P)}. (3.42b) 
with the anomalous clastic constants 

C^ii) = C2r^Z(^)-3/2 [t{e)/t]-' [aii)/af^ , (3.43a) 

Csii) - c^r^zii)-' im/tr' He)/^r'^' , (3.43b) 

Ci{i)^C3ie). (3.43c) 

To obtain the equation for C2{i) we used the fact that 
it must be proportional to C2 and thus to a^^^. Note 
that a affects the leading behavior of C2 despite flowing 
to zero. In other words: cr is a dangerous irrelevant vari- 
able. Similar arguments for C3 imply that C3{£) must be 
proportional to p which reaches a non-zero fixed point 
value. There are, therefore, no contributions to the scal- 
ing of C3 from dangerous irrelevant variables. 



Combining the formal scaling form of Eq. H3.35|l with 
the power laws of Eq. H3.45|l for Z{i), lu{£)/lu and d-{£)/a 
and the crossover length scales as given in Eqs. (|3.44|) 
and (|3.37b|) . we obtain a complete picture of the scaling 
behavior of the displacement vertex functions in d < 3. 
In particular we obtain Eqs. H2.2|) for the 2-point func- 

featuring the 



tions r 



(0,2) 
dd 



scaling exponents 



riy^ and rif) with 



?7i^ =£ + 7* + C72 = 4£/7, 
?7c = -e-27*-~C72 = £/7, 
m = -£-3Y/2 + w/2 = e/7 
= 2-C72 = 2-2e/7 = (4 - 



(3.46a) 
(3.46b) 
= f]c , (3.46c) 
??K)/2.(3.46d) 



Thus, even though there are three independent expo- 
nents 7*, and w the scaling behavior of the 2-point 
functions is determined, at least to first order in e, by 
only two exponents, say rjK and rjc- Upon inserting the 
power laws (|3.45f) into Eqs. H3.43|l and (|3.40a|) we obtain 
the scaling expressions for C2, C3, C4 and K stated in 
Eqs. H2.5() and H2.6|l the review section. Note that K 
diverges at long length scales whereas C2, C3, and C4 
vanish in this limit. 



5. Behavior of the elastic constants in d < 3 

For e > we can assign a physical meaning to the 
hitherto arbitrary length scale L± via the definition of 
our dimensionless temperature t, viz. 



Li = 



(3.44) 



Of course, this choice also affects the length scale Ld, 
cf. Eq. (|3.37bp . The length scales L± and Ld mark the 
borderline between harmonic and critical behavior. 

Below 3 dimensions the solutions to the characteris- 
tics are for ^ ^ 1 governed by the IR-stable fixed point. 
Readily, one finds the power laws 

Z{£) = , (3.45a) 

uii) = uif , (3.45b) 

where 

7* = 7(i*,f^*,P*)-~5£/7, (3.45c) 

C* - C{t\(J*,P*) ='^ell . (3.45d) 

From our discussion above it is clear that we cannot sim- 
ply set a{€) equal to its fixed point value cr* = because 
(7 is a dangerous irrelevant variable that effects the scal- 
ing behavior of C2 at leading order. The vanishing of a 
is described by 



(3.45e) 



with a Wegner exponent w — ejl corresponding to 
the smallest eigenvalue of the Hessian of the fiow near 



6. Logarithmic behavior in d = 3 

Since e vanishes in d = 3 the solutions to the character- 
istic equations are no longer of power law type. The fiow 
of the temperature, for example, is described at leading 
order, i.e., for a{£) = a* and p{£) = p*, by 



7 



(3.47) 



32 TT 

This differential equation is readily solved with the result 

-1 



m/t 



1 



7t 
327 



In(^) 



(3.48) 



Similarly, we find Z{i) ~ [t{i)/t]-^, u;{e) - [t(^)A]^ 
and a{£) ~ [i{£)/t]^ at leading order. Inserting these 
logarithmic solutions into Eqs. H3.40a|l and l|3.43|l we find 
that 



C3(qj.) 



C2 



1 



1 



It 
32^ 

It 



ln(i±|qil) 



32 TT 

CMi_) - C2(q 



ln(i±|qxl) 



4/7 



-1/7 



(3.49a) 

, (3.49b) 
(3.49c) 



Of course, Ci is normal as it was in d < 3. Once more, K 
diverges at long length scales and C2, C3, and C4 vanish 
in this limit. 

Note that Eqs. (|3.49|) imply the existence of a nonlinear 
crossover length scale, viz. 



Cj. = L± cxp 



32 TT 

TT 



= cxp (^^^j. (3.50) 
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For C^lqj^l « 1 the anomalous elastic constants are ap- 
proximately harmonic whereas one has clearly anomalous 
behavior for C_L|q_L| <C 1. 



7. Cs as a relevant perturbation - semi-soft elasticity 

Up to this point we have excluded a term proportional 
to C5, cf. Eqs. (|2.1|) and (|3.3|) . from the elastic energy 
of CSEs because such a term destroys the soft elasticity. 
Now we establish contact to more conventional uniaxial 
elastomers by incorporating a small but nonvanishing C5. 
As a consequence, we find semi-soft behavior. 

Technically, we treat C5 as a perturbation to the CSE 
model. It turns out that this perturbation is relevant in 
the sense of the RG. This situation is analogous to the 
(j>* model where a deviation from the critical temperature 
represents a relevant perturbation. Our central task we 
will be to determine the scaling exponent that governs 
the departure of C5 from zero. 

Our analysis here is based on the full Hamilto- 
nian (|2.1I) . Carrying out the /i-rescaling it is straight- 
forward to see that not the entire strain Uad is relevant 
and that it is sufficient to keep 



Uad = ^daUd ■ 



(3.51) 



Applying the rescaling that led us from Eq. (|3.4|) to (|3.t)|) 
we obtain 



He 



d''^x± / dxddaUddaUd, (3.52) 



with H as stated in Eq. H3.6(l and where e = C5/{4:K). 
Note that e ^ fj?, i.e., the naive dimension of e is 2 and 
hence e is clearly relevant. 

The Gaussian part of Tie has an extra term compared 
to Ti and hence the Gaussian propagator of He is different 
from that of H. Here, any in Eqs. (|3.12(l has to be 
replaced by q^ + eq^. The non-Gaussian terms of He 
and Ti are identical and hence we still have the 4 vertices 
stated in Eqs. (EUSJ- 

To investigate the departure of e from zero we expand 
the propagator to linear order in e. Then this expanded 
propagator is used in our diagrammatic calculation. Of 
course, at zeroth order in e we retrieve our vertex func- 
tions H3.14|) . The first order in the expansion leads to 

an extra divergent term in E^j'^''^'' that is proportional to 
e q^^ . To remove the extra divergence we introduce an 
additional rcnormalization factor via setting 

e ^ e — Z^^ZeC . (3.53) 

^From the diagrams depicted in Fig. ^ we extract that 

6 + 5p + p2 



16ne^/2Tp^/2 + p- a 



(3.54) 



The RGE for the vertex functions expanded to linear 
order in e reads 



D„ 



M+|)7 



X r(*^'^) i{q±,qd} ; LU, t, a, p, e,^^) = (3.55) 
with Df^ as stated in Eq. H3.22|l and 

V' = /i9^1ne|o = 7t - 7e ■ (3.56) 
Setting up a characteristic for e, 

^^lne-(£) = i^m^m.m) , e-(l) = e, (3.57) 
we find that this coupling flows in d < 3 as 



(3.58) 



where ■0* = il^{t* , cr* , p*) — — e/7. Now, the solution to 
the RGE H3.55|l in conjunction with dimensional analysis 

tells us that the scaling behavior of F^;'^''^'' expanded to 
linear order in e is given by 



r 



i"/^ {q±,qd;u;,t,a,p,e,fi) = i^le)^-' 2(1)-^ (3.59) 



Next we switch back to physical variables. By performing 
much the same steps as in Sees, pil B 5ll and B 4p we 
obtain for c? < 3 



^dd i^±^<ld;e) 



K 



(3.60) 



X <i> 



(0,2) f L ±q± Ldqd, L\e 



dd 



where we have dropped several arguments for notational 
simplicity and where 



1/5 = = l/2-£/28. 



(3.61) 



Thus, C5 plays the same role in this problem as tem- 
perature plays in a traditional thermal phase transition. 
As in the thermal case, it is useful to introduce a corre- 
lation length 



^5 = L^ieLlY 



(3.62) 



There are two interesting limits we can now consider: 

^5|q_L| < 1 and <^5|q_L| > 1- In the first case, F^^'^^ 
must be proportional to q^ when = and to when 
qj^ ~ 0. Thus, we obtain by choosing I — {£,5/L±)^'^^ 
from Eq. ^W^: 



-n(0.2) / X 

^dd (q-L,9<i;e) 
T I cJd 



If KLl^-^-^'-'^'^^ql^Crdl if qd = 

if q 



(3.63) 
) 
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where 

75 = 1^5(2 -7?A') = l-5£/14. (3.64) 

The second hmit ^slqj.! S> 1 corresponds C5 0, and 
we must obtain the same scahng forms as we obtained 
for C5 = with a correction term that vanishes with 

As we stated earher, we have not inchided nonhnear 
terms in our model Hamiltonian [Eq. 13. 4|) ] that are 
needed to stabilize the system when C5 < 0. In particu- 
lar, we have not included terms {uadUad)'^ ~ {daUddaUdY , 

UaaUbdUbd, UabUadUbd, ula^bdUbd, ul^UcdUcd, and 

UaaUbc'^bd'^cdi ^11 of which havc the same naive scaling 
dimension as the harmonic terms in nonlinear strain 
that we retained in Eq. (|3.4() . Our expectation is that 
the general theory in which these terms are included 
will have the same general form as the present theory 
with, however, a different stable fixed point with in 
particular a nonzero value of the coefficient of [uadUadY ■ 
At such a fixed point, (daUd) will develop a nonzero 
value at negative C5 that scales as {—C^)^ where 
/3 = (1 + jj*)'^5- This result can be obtained by observ- 
ing that Eq. for T^'''^) implies d^Gadixd, a;a, C5) = 
l2+i'/2gaQ^^(^l2+C/2,.^j^^j-^.C5), where d', 
is a derivative with respect to Ixa, and that 

Gdd{xd,Xa,C^) {daUd){daUd) aS X-+ CX3. 



IV. NEMATIC ELASTOMERS 

Usual NEs are either crosslinked in the isotropic or in 
the nematic phase. If synthesized in the isotropic phase, 
the uniaxial anisotropy arises via a spontaneous symme- 
try breaking at the isotropic to nematic transition and 
is associated with soft elasticity. Crosslinking in the ne- 
matic phase, on the other hand, permanently imprints 
the uniaxial anisotropy into the material and leads to 
semi-soft behavior. We will start by studying the soft 
case. Further below, we will include the effects of an im- 
printed uniaxial anisotropy to investigate the semi-soft 
case. 



A. The model 

Since the spontaneous symmetry axis of soft NEs can 
point in any direction, their elastic energy has to be ro- 
tationally invariant not only in target space but also in 
reference space. Both invariances are taken into account 
by writing the stretching energy as 




+ A2 tr u. tr + ^3 tr + Si ( tr ^) + B2 ( tr 

B3 ( tr u) ^ tr + B4 tr u tr + B5 tr I , (4.1) 



with the first two expansion coefficients being the usual 
Lame coefficients. Of course, terms of higher than fourth 
order are allowed by the symmetries of the system. How- 
ever, these higher-order terms turn out to be irrelevant 
in the RG sense and are hence neglected. 

Suppose that the spontaneously uniaxially ordered 
elastomer is described in equilibrium by an equilibrium 
strain tensor u^. Without loss of generality we may as- 
sume that the anisotropy axis lies in the direction and 
that is a diagonal matrix with the diagonal elements 
uoaa — uo_L and UQdd = uo\\- To describe deviations from 
the equilibrium configuration, wc introduce the relative 
strain 

w_ = u — u^ (4.2) 

and to expand TH. in terms of w_. By dropping terms that 
depend only on we find 

Hst^ j d'^^xx_ j dxd{aiWdd + a2Waa + bi w}^ 

+ h WddWaa + h + &4 + 65 + Ci Wddwl^ 

+ C2 Waawld + ^3 WabWadWbd + di W^^m} 7 (4.3) 

where we have discarded terms that turn out to be irrel- 
evant. The new coefficients ai, 02, &i, and so on, depend 
on the old coefficients A, fi, and so forth, as well as on 
Mo II and uo±. 

By virtue of the rotational invariance in reference space 
there exists a set Ward identities relating the vertex func- 
tions implicit in H4.3|l . We derive these identities in Ap- 
pendix^ At zero-loop or mean-field level these Ward 
identities correspond to relations between the elastic con- 
stants in (|4.3|) . 

ai — a2 — s 65 = , (4.4a) 
62 - 263 - s C2 = , 62 + 65 - 26i + s ci = , (4.4b) 
65 - 2&4 - s C3 = , ci - C2 - C3 - 2s di = , (4.4c) 

where s is an abbreviation for Ug|| — Uo_l- Since w de- 
scribes deviations from the equilibrium its thermal av- 
erage (w) has to vanish. At zero temperature, where the 
mean- field approximation becomes exact, this means that 
the coefficients of the linear terms in (|4.3f) must be zero. 
Equation H4.4a|l then leads to the observation that 65 = 
for uo|| uo±. At finite temperatures thermal fluctua- 
tions become important and loop corrections renormalize 
the elastic constants including ai, 02 and 65. For (w) to 
vanish the renormalized versions of oi and a2 have to sat- 
isfy equations of state to which ai — and a2 = are the 
mean-field approximations. In the following we will as- 
sume that we have chosen ai and a2 appropriately so that 
their respective equations of state are satisfied. In other 
words: we assume that we expand about the true equilib- 
rium state. Then, the Ward identity (|Blip generalizing 
(|4.4a(l guarantees that the renormalized 65 vanishes for 
""oil 7^ ^^O-L- The vanishing of the elastic constant 65 is 
the origin of the softness of NEs as can easily be seen by 
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rewriting the nonvanishing terms of the Hamiltonian at 
leading order in small deformations in Fourier space. 

In what follows we assume that ai = a2 = 65 = 0. 
Exploiting the relations (|4.4|l we rewrite the stretching 
energy as 

Hst = J d'^^X±_ j dXd {bi + 62 VddVaa + h vl^ 

+ b4vl,}, (4.5) 



where we have introduced the non-standard strains 



Vab 



Wab - 
Wdd 



WadWbd , 
^WadWad ■ 



(4.6a) 
(4.6b) 



Next we cast Vab and Vdd into a more familiar form. Re- 
call that the relative strain w_ depends upon the displace- 
ment u relative to the original isotropic state measured 
in the original reference space coordinate x. It is more 
convenient, however, to work with a relative strain v! 
that depends on the displacement u' relative to the equi- 
librium state of the NE measured in the coordinate x' of 
the corresponding uniaxial reference state. The tensors 
w and u' are related via (see, e.g., Ref. [Tlj |) 



(4.7) 



where ^ is the so-called Cauchy deformation tensor of 
the uniaxial equilibrium state. A^^ and provide equiv- 
alent descriptions of this state and they are related via 

io-^(4"4-i), (4.8) 

where 1. denotes the d x d unit matrix. Substituting the 
relation H4.7|l into Eqs. H4.6|l we obtain 



Vab 



Vdd = '"Aoj^ 



d'bU'a - 



1 1 



^a"'d"a'^d 



2 r - 1 



, (4.9a) 
(4.9b) 



where we have used Eqs. (|3.5|l and (|3.51|) to express u'. 
r = Ap||/Aq^ is the usual anisotropy ratio that charac- 
terizes the anisotropy of the uniaxial equilibrium state. 
In the steps leading to Eqs. H4.9|l we have exploited that 
s = Aoj^(r - l)/2. A glance at Eqs. H4.9() shows that 
we can write w in a simpler and more traditional form 
by rescaling x'^ Xa, x'^ ~> \/r ~ 1 Xd, u'^ — > Ua and 
u'j^ —> -^/r — 1 Ud- Incorporating bending, we eventually 
arrive at the Hamiltonian 

2 

+ 2 C2 VddVaa + C3 vl, + 2 C4 «^ J , (4.10) 
with the non-standard strains 

Vab = \{daUb + dbUa - daUddbUd) , (4.11a) 



n 



d'^^x^ / dxd{C^vl^ + K{Vludy 



Vdd 



ddUd + \daUddbUd ■ 



Here, we exclusively retained the only bending term 
that is relevant in the RG sense. The new elas- 
tic constants C are proportional to the h in 
^■.Ci = 2 6iA4^r2VF^, C2 = 62A4^r^/7^, C3 - 
2 bsAfj^ and C4 = hik'^^^/T^. 

A comment regarding the rescaling of Xd is called for. 
Our rescalings x'^ — > \/r — \xd and u'^ — > \/r — \ud only 
make sense if r > 1. We have assumed that the ne- 
matic phase is characterized by r > 1. Nematic phases 
with r < 1 are also possible. In this case, our rescalings 
would be different. At a second-order phase transition, 
r — 1 will tend to zero. The transition from the isotropic 
to the nematic elastomer phase is first order in all di- 
mensions above 2. The order of the transition in two 
dimensions has not yet been established. A second-order 
transition would have unusual properties since the coef- 
ficients of the nonlinear contributions to the nonlinear 
strain of Eqs. (|4.9I) diverge as r ^ 0. 

Instead of using the Hamiltonian 1)3. 4|l or our RG anal- 
ysis we find it convenient to reduce the number of con- 
stants featured in the statistical weight exp(— 7i/T) by 
^/K^4T, Xd 



rescaling T 
^jK/dud and Ua 



^JCi/Kxd, Ud 
{K/C4)ua- This gives us finally 



(4.11b) 



T " Jf j '^''^^^ j dxd{Lovli+{Vludf 

+ 2gvddVaa + fvl, + 2vl^,], (4.12) 

where 

c^ = Ci/C4, .g = C2/C4, / = C3/C4. (4.13) 

Like the parametrization of Eq. H3.6|l of the Hamiltonian 
for CSEs, the parametrization of Eq. (|4.12|l is not appro- 
priate for taking the C4 ^ 00 limit to obtain a smectic-yl 
Hamiltonian. As was the case for CSEs, we are interested 
in properties unique to NEs, and we will not consider the 
Grinstein-Pelcovits limit of our model. 

Formally, the Hamiltonians (|3.6|l and (|4.12|) look very 
similar. One has to bear in mind, however, that the 
strains u and w_ are different. In fact, the scaling 
symmetries of (|3.6|l and (|4.12|) are quite different ex- 
cept for the ^ invariance. The naive dimensions of 
the fields, the temperature and the coupling constants 
are the same for the two models. In particular, NEs 
and CSEs have a mutual upper critical dimension dc = 
3. Though both systems possess of a mixing invari- 
ance the specific forms of these invariances are distinct: 
the NE Hamiltonian is invariant under the transforma- 
tion Ua{Xc,Xd) Ua{xc - O^XdyXd) + 9aUd{xc,Xd) and 
Ud{xc., Xd) Ud{xc ~ OcXd, Xd) + OaXa provided that the 
0's are small. This mixing invariance is reminiscent of 
the original reference space rotation symmetry. In our 
current approach to NEs, the /3-invariance of CSEs has 
counterpart. However, in Appendix IdI we present an al- 
ternative formulation that features a /3-invariance at the 
cost of having an extra scaling parameter. Due to the 
different forms of the strains, the remaining symmetries 
stated at the end of Sec. Illl^ have no analog in NEs. 
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B. Renormalization group analysis 



2. Renormalization 



The diagrammatic perturbation expansions for CSEs 
and NEs are similar. Instead of repeating the details 
we will high-light the differences. In contrast, the RG 
behavior of CSEs and NEs is quite different due to the 
varying scaling symmetries. 



1. Diagrammatic expansion 

It is not difficult to see that the Gaussian propagator 
for NEs coincides with that for CSEs, see Eqs. I|3.12|l . 
The differences between u and v, however, lead to differ- 
ent vertices. For NEs we have 



. ^ -g (1) (2) (3) 
'^^2Tld lb % ^ 

■ 9- f (1) (2) (3) 

■ 1 (2) (1) (3) 

- 2ff + / + 2 (1) (2) (3) (4) 
Ha % % 



(4.14a) 
(4.14b) 
(4.14c) 
(4.14d) 



Of course the sum of the momenta has to vanish at each 
vertex. Note that the vertices (|3.13|) and (|4.14|) are of the 
same structure. Merely the coupling constants appear in 
different combination. Hence, the Feynman diagrams for 
both models have the same topology, or in other words, 
the 2-leg diagrams for both models can be drawn as in 
Figs, n to 131 Moreover, the same type of vertex func- 
tions are superficially divergent. By virtue of the NE 
mixing invariance, there exists a set of Ward identities 
relating the NE vertex functions. These identities are 
derived and stated in Appendix |21 They guarantee that 
we merely have to calculate the 2-point functions. Ap- 
pendix ^describes details if the reader makes the appro- 
priate changes in the vertices. Our results for the 2-point 
functions read 



(4.15a) 



16 7r£ V-5' + (2 + /)^ 
_ [Ag + g^-2{(S + uj)~f{l2 + uj)\ gj ^^^/^ 

32 7rex/2T7V-5' + (2 + /)c^ 
ri'/^ (4.15b) 
_ V2 / (1 -I- / - g) [g ~ u) qj. Qa ^_^/4 
16ney/~g^ + {2 + f)u; 

r['f ^T-'[{f + l)q,q, + SabCil] (4.15c) 



V2T7[2{1 + f^grqaqb + Sab <ll] 

T 



/4 



We eliminate the e poles from the NE vertex functions 
by employing the renormalization scheme 



Xd - 


^ id 


1 /n 

= Z-^I^Xd. 


(4.16a) 


Ud - 


* Ud 


^Z^/^Ud. 


(4.16b) 


Ua - 


■> Ua 


= ZUa , 


(4.16c) 


T ^ 


f = 




(4.16d) 


UJ — > 


ill = 


- Z Zrp Z ^ijj , 


(4.16e) 


9 


9 = 


Z ^Zrj?^Zgg^ 


(4.16f) 




/ = 


Z ^Zrj^^Zfj . 


(4.16g) 


tions 


are 


devised so that the strains Vab 



and Vdd as well as our Hamiltonian remain invariant in 
form, 



T 2~ 



^ j d'^^xi_ J dxd[z^ujvla + Zt {V\udf 

+ 2ZggVddUaa + Zffvl^ +2ZTZVabVab^ , (4-17) 



The scheme (I4.16f) follows closely the approach devel- 
oped by Grinstein and Pelcovitz Of course, other 
re-parameteriztions are conceivable. In Appendix ^ we 
present an alternative formulation with a different renor- 
malization scheme in which neither the elastic displace- 
ment nor Xd is renormalized. 

In our current formulation, there remanis no scaling 
invariance of the Hamiltonian that can be exploited to 
further reduce the number of coupling constants. Hence, 
the renormalization factors are functions of the original 
dimensionless parameters t, uj, g and / rather than of a 
reduced number of effective couplings. The structure of 
the NE renormalization factors is 



Z,„{t,uj,gJ) = l+Y. 



X(")(i,c^,g,/) 



(4.18) 



with X*^"')(t, w, g, /) being a power series in the effective 
temperature t beginning with the power i™. To 1-loop 
order we find from (|4.15|l via minimal subtraction 



Z =l + t 



4,g + ff^-2(7-Fc^)-/(13 + cj) 
32 Tre V2T7 ^-9^ + (2 + /) w 



\-t 



Ag + g"^ -2{% + Lo)- f [12 + uj) 



Z^ = l + t 



Zn = l+t 



Zf = l + t 



32 7reV2T7V-.9' + (2 + /)^ 

V2Tfig-Lof 
mTTSujy/-g^ + (2 + /)tJ ' 
V2T7(l + /~g) jg-Lo) 
16 7regV-.9' + (2 + /)^ ' 
^/2T7[l + 4(/-g) + 2(/-ff)2 
32nsf^-g^ + i2 + f)u; 



32Tre^-g^ + {2 + f)i 



(4.19a) 
(4.19b) 
(4.19c) 
(4.19d) 

(4.19e) 
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3. Scaling 1: RG equation and its solution 

The RGE for the NE vertex functions follows as ususal 
from the fact that the unrenormalized theory has to be 
independent of the arbitrary length scale introduced 
by renormalization. Instead of working with the original 
parameters lj, g and / we prefer to switch to 



K = g/io = C2/C1 , 
(j=1/lo = C4/C1 , 



p = //c^ = C3/C1 , (4.20a) 
(4.20b) 



This step turns out to be helpful in studying the RG flow 
because some of the original parameters tend to flow to 
infinity. We will see shortly that, on the other hand, k, 
p and cr flow to finite values. Our RGE reads 



N 



7 r(*^^^)({q^,g4;t,«,p,(7,/i) = 

(4.21) 

with the RG differential operator 

(4.22) 

The Wilson-/3-functions, from which we determine the 
fixed points of the RG flow, are given in terms of the 
Wilson-7- functions 



7... 



^ 1 z 
fi— InZ.. 

0/1 



by 



(3t = i (-£ + 7T - 7/2) : 

Pk. = K iluj - Ig) , 

Pp = PiliJ- If) , 

Pa = O- {juj - 7T - 7) • 



(4.23) 



(4.24a) 
(4.24b) 
(4.24c) 
(4.24d) 



The Wilson-7-functions are readily extracted from the 
renormalization factors with the result 



j,„^~tdtX(]\t,u;,gJ). 



(4.25) 



Switching to the parameters defined in Eqs. H4.20|l we 
obtain 



Pt - -te 



+ 12Kcr - 2(t(3 + 19a) 



(4.26a) 
p(3 + 37a) 



64 TT cr ^JYVpja p + 2a 



P. 
Pp 



v/2 + p/a (k - 1) (p + cr - K^) 
16 7r + 2cr - k2 



^2 + P/^ 



32 TT + 2cr - k2 
X {2k^{1- p) + 2p^ -Ana + a"^ 



p{Aa^ 



Pa = t 



^2 + pjaa (-2 + 4k - 2k^ + a) 

32 TT + 2(7 - k2 



(4.26b) 

(4.26c) 
2)} 

(4.26d) 



For solving the RGE we introduce the characteristics 



dl 



p, p(l) = p, 



(4.27a) 



i^^\nz{t) = im,m,mrpm , ^(i)-i 



(4.27b) 
t(l)=t, 

(4.27c) 

fii(l) = K , 

(4.27d) 

pW = p, 

(4.27e) 
a{l)=a. 

(4.27f) 

In contrast to the CSE model we have to look for fixed 
points in a 4 dimensional parameter space. We find that 
the quadruple of coupling constants {t{tj, K{i), p{l),a{£)) 
flows to the IR stable fixed point 



eg-^m = ptm,m,m,m) 



i-m = PAm,m,m,m) 



i-m = Ppim,m,m,m) 



i-a{e) = pAiii),m,m,m) 



(4.28) 



This fixed point is characterized by C2/C1 = k* = 1 and 
C3/C1 = p* = 1. It turns out that the leading scaling 
behavior of physical quantities depends not only on the 
fixed point but also on the approach to the fixed point 
described by the dangerous irrelevant variable a. Paths 
to the fixed point decay quickly to the line described by 



n{l)~K* = (a/2 + l/4)o-(f) 



(4.29a) 



and 



p{t) - p* ^ aa{t) (4.29b) 

for small where a is an arbitrary constant. In addition 
to the stable fixed point, there is the unstable Gaussian 
fixed point t* = Q and there are two unstable fixed lines 
which can be parameterized as (i(cr), K(cr), p(cr), a) with 



t{a) 



64 7r£ \/2~~272aTia 



3136 - 3320cr + 1521cr2 
X |56V2 + 124^/^+39\/2o•| 

k(ct) = 1 - ^/aJ2, 
p(cr) ^ l-\/2a-cr/2 



(4.30a) 

(4.30b) 
(4.30c) 



and 



t{a) 



64 7re\/2 + 2\/2CT + 3cr 



3136 - 3320cr + 1521cr2 
X |56\/2- 124V^+39\/2(t| 

k(ct) = 1 + v/a72, 
p(cr) = l + V2^ -a/2. 



(4.31a) 

(4.31b) 
(4.31c) 
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With help of the characteristics a formal solution to 
the RGE is easily obtained, 

r('^'^) ({qx, 9.} ; t, «, p, a, m) = (4.32) 
X r(^^'^) ({q^, t-(^), p{e),a{i), M^) ■ 

Dimensional analysis gives 

r(^^'^) ({q^, g.} ; t, n, p, a, p) = ^-(^+i)+^m+(.+i)a^ 
xr(*^'^)({^,^};t,«,P,^,l). (4.33) 

Equation 1)4. 32|) in conjunction with (|4.33|) finally reveals 
the scaling behavior of the vertex functions: 



(4.34) 



4-. Scaling 2: physical quantities 

Now we switch back from the convenient scaled vari- 
ables that we used in our calculations to the original 
variables in which we formulated our Hamiltonian H4.10|l . 
Once more, we mark rescaled variables as well as vertex 
functions of the rescaled variables with hats. Recalling 
our manipulations leading to Eq. (|4.12(l we write 



T 




T, t 



C4 C*4 
Id = -f7<ld , Ua = -frUa , Ud 



(4.35a) 
'^Ud. (4.35b) 



The relation between the physical vertex functions and 
the vertex functions in the scaled variables is given by 

r(*^'^)({q^,g4;i,K,a,p,M) 

^ ^Af+37V/2-7/2^3/2-JV/2 

xf(*^,^)({q^,g4;£,/t,fT,p,M) . (4.36) 
Equation (|4.36() in conjunction with Eq. (|4.34() gives 

T'^^''^H{q^,qd};t,K,a,p,p) (4.37) 

= }_Ty'M+3N/2-2(^l-N/2j^2{d-l)~dM-{d+l)N 

4 _L 

X i-(d+l)+dM+(d+l)N ^t{l)/t]-^ Z(^)-(^'^+l)/2 



X ^iM,N) 



ij.qj_ Ldqd 



£ ' Z(£)l/2£2 



^tii),m,ai£),p{e)a 



where we introduced the susceptibilities 



(4.38) 



and where we switched from p to the length scales 

= p-^ , (4.39a) 
Ld-Ml^-'-MLl. (4.39b) 



The behavior of the elastic moduli is now easily ex- 
tracted from the 2-point vertex functions. Upon special- 
izing Eq. to iV = 2 and Af = etc. we obtain 



(4.40a) 



C,{£) = C,£-^Z{£)-y^ im/tr' [a{£)/'Tr' , (4.40b) 



Ci{£) = C4£-'z{£)-^^^ [t{e)/tY 

C2{£) - C3{£)^Ci{£). 



(4.40c) 
(4.40d) 



For Ci, C2 and C3 we can draw our final conclusions 
without solving the characteristics. Taking the derivative 
with respect to the flow parameter we find, for example, 



£§-^CAi) 



_3 

2'^ t a 



(4.41) 



Upon expressing in terms of cr, p, and k and by taking 
into account Eqs. H4.29|l we see that 



7c. 



(4.42) 



in the vicinity of the IR-stable fixed point. Hence, 
vanishes in the scaling limit. Consequently, Ci, C2 and 
C3 are normal. 



5. Behavior of the elastic constants for d < 3 

For d < 3 we can assign physical content to the length 
scale L_L by exploiting the definition of the rescaled di- 
mensionless temperature t. This provides us with 



l/e 



(4.43) 



as the transversal length scale associated with the on- 
set of scaling behavior. The corresponding longitudinal 
length scale can be inferred from Eq. ()4.39b|) . 

As solutions to the characteristics we obtain the power 
laws 

Z{£) ~ r* , (4.44a) 
d{£) - ar , (4.44b) 



where 



and 



7* 7(r , K*, p*, cr*) = -42£/59 (4.44c) 



u; = 4£/59 (4.44d) 
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is a Wegner exponent corresponding to the smallest 
eigenvalue of the Hessian at the stable fixed point. 

Equation (|4.37(l in conjunction with the power-law be- 
havior of Z{£) and ^{i)/(7 along with Eqs. H4.43|l and 
llOObll for the crossover-lengths L± and Ld provide us 
with a complete scaling picture of our vertex functions in 
d < 3. This picture is summarized in Eqs. H2.13|) with 

r]K = s + J* /2^38e/59, (4.45a) 
r]c = s + 3-/*/2^Ae/59 = w, (4.45b) 
= 2-1-772 = 2- 2l£/59. (4.45c) 

Finally, the power-laws (|4.44|l along with Eqs. H4.4()|l re- 
sult in the scaling forms for the elastic constants sum- 
marized in Eqs. H2.16|l and H2.17() . Note that K diverges 
at long length scales. C4, on the other hand, vanishes in 
this regime. 



6. Logarithmic behavior in d — 3 

In 3 dimensions the solutions to the characteristics that 
we need to determine the behavior of the elastic moduli 
are given at leading order by 



m/t 



7 Vet 

647r 



In(^) 



(4.46) 



and Z{£) ^ [t{£)/t] 59 . Inserting these logarithmic solu- 
tions into Eqs. (|4.40|l we obtain 



C4(q^) ^ Ci 



7^6 1 
647r 

7 Vet 

~ 647r 



ln(L±|qil) 
ln(i^|q^| 



38/59 

, (4.47a) 

-4/59 

. (4.47b) 



Using (t(£) ~ [f(£)/i]59 one can check explicitly that Ci, 
C2 , and C3 are normal. K and C4 diverge and vanish, re- 
spectively, for ij^lq^l 0. Note that Eqs. (|4.47|l imply 
the existence of the crossover-length scale 



C_L = Li_ exp 



/ 647r 



Li_ exp I ^ I . (4.48) 



C4/C1 = 0. Remarkably, not only the stable fixed point 
but also the approach to it contains information on re- 
lations between the elastic moduli. From Eqs. (|4.29|) it 
follows that 2k — p— l = cr/2 for small flow parameter £. 
Switching back to the original elastic constants we find 
the Poisson ratio 



2 C2 — C3 — Ci 



(4.49) 



Note that Eq. H4.49I) implies the Poisson ratio h/hl — 8 
found by Xing and Radzihovsky p^ . where fj.^ and /i 
are longitudinal and transverse shear moduli, respec- 
tively [321 ■ 



8. Semi-soft elasticity 

So far we considered soft elastic NEs synthesized by 
crosslinking in the isotropic phase. If a NE is crosslinked 
in the nematic phase, a memory of the anisotropy at the 
time of crosslinking is locked in, or in other words, the ro- 
tational symmetry in the reference space is broken. The 
simplest way of modeling this symmetry breaking is by 
introducing an aligning external field, or more precisely, 
an aligning external stress. In this spirit we supplement 
the elastic energy (|4.1() with the term stated in Eq. H2.11(l . 
Next, we switch to the relative strain 14.2|l and exploit a 
Ward identity that corresponds at 0-loop order to 



h — 2 sbc. 



(4.50) 



By virtue of this Ward identity, which is derived Ap- 
pendix we can study the effects of the external align- 
ing stress via studying the RG behavior of 65 featured in 
Eq. g2J|. 

Before embarking on a perturbation calculation we re- 
call the rescalings that led from Eq. ()4.3|l to our final 
NE Hamiltonian H4.12|l . By applying these rescalings to 
Eq. (|4.3|) with now — h/s instead of 65 = 0, by incor- 
porating bending and by dropping irrelevant terms we 
derive the model Hamiltonian 



He 

T 



n 

T 



e 
2f 



Xl_ / dXddaUddaUd, (4.51) 



with H/T given by Eq. (|4.12|) . The coupling constant e 
is defined as 



7. Poisson ratios 

Since the momentum dependence of the anomalous 
elastic moduli is logarithmic in 3 dimensions it will be 
difficult to observe the anomalous elasticity of NEs upon 
measuring these moduli directly. However, our analy- 
sis reveals the existence of several Poisson ratios of the 
elastic moduli that should be conveniently accessible by 
experiments. 

The IR stable fixed point values of k, p and a directly 
imply the Poisson ratios C2/C1 = 1, C3/C1 = 1 and 



4i^ (1 - 2s) 



(4.52) 



Being interested in semi-soft behavior we assume that h 
is small. Hence, we treat e as a relevant perturbation. 

Our diagrammatic calculation proceeds in much the 
same way as the one outlined in Sec. lIIIB7l The only dif- 
ference is that we here have to replace the vertices (I3.13f) 
by the NE vertices H4.14|l . Expansion of the diagrams to 

first order in e leads to an e pole in F^^'^^ that is propor- 
tional to eq\. We remove the divergence by letting 



(4.53) 
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with 



CONCLUDING REMARKS 



x/2+7(3 + /-2g + cj) 
\QTie^-g^ + {2 + f)uj 
The appropriate RGE here reads 



(4.54) 



D,. 



e ip de 



N 



X r(*^^^) ({q^, g4 ; t, K, p, a, e, Ai) - (4.55) 

where it is understood the vertex functions are expanded 
to first order in e. The RG differential operator is 13^ is 
stated in Eq. H4.22|l and 



-0 = fidf^\ne\o = 7t 
In d < 3 the characteristic for e, 
d 



(4.56) 



^P{m,m,m,m) , e(l) = e,(4.57) 

has the fixed point solution 

e(^)-e^'^*, (4.58) 

with =ilj{t*,K*,p*,a*) = 18e/59. Supplementing the 
solution to the RGE H4.55() with a dimensional analysis 
and switching to original variables yields for the physical 
vertex function ^^jif^ the scaling form 



K 



X $ 



(0,2) 
dd 



T 



(4.59) 



where we have simplified the notation by dropping several 
arguments and where 



= (2-1/'*)"^ = l/2 + 9e/108. 



(4.60) 



We emphasize the different roles played by C5 and h in 
the CSE model and in NEs, respectively. C5 corresponds 
to a temperature whereas h corresponds to an external 
magnetic field. Nevertheless, it is useful to introduce a 
length scale 



at this point. In the limit .^hlqj 
choosing £ = {£_h/ L^)-"'- that 



' h-"'^ (4.61) 
<C 1 we obtain upon 



-n(0,2) / X 



1 



s/i q_L 



if 
if 



(4.62) 

qd^Q 

q^ = 



where 



7^ = z/;,(2-77x) = l-10e/59. 



(4.63) 

In the limit ^/i|qi| 3> 1, corresponding to the limit 
/i — > 0, we retrieve essentially the scaling form that we 
had for h — 0. Now, however, there is an additional cor- 
rection term that vanishes as {^h\n±\)~^^'^'^ ■ Physically, 
this correction term modifies the behavior from soft to 
semi-soft. For a review of the complete scaling results, 
see Eqs. ^J^, and (|TT7|) . 



In this paper, we have explored the anomalous elas- 
ticity of two models for soft uniaxial elastomers, both 
of which are characterized by the vanishing of the shear 
modulus C5 for shears in planes containing the anisotropy 
axis. The first model, which we refer to as a critically soft 
elastomer (CSE), describes a uniaxial system at a simpli- 
fied critical point separating a true uniaxial elastic phase 
characterized by five elastic moduli and a lower symme- 
try phase produced by shearing the uniaxial solid. The 
second model describes nematic elastomers (NEs) formed 
via spontaneous symmetry breaking from an anisotropic 
state. In the CSE model, a bending modulus diverges, 
and three of the four clastic moduli vanish as power-laws 
in wavenumber at long wavelength below spatial dimen- 
sion c? = 3 and logarithmically at d = 3. In the NE 
model, the relevant bending modulus diverges, but only 
one elastic constant exhibits singular behavior. In both 
models, we studied the effect of turning on couplings that 
take them from soft to true uniaxial elastomers. These 
fields act like temperature or external fields at a ther- 
mal critical point, respectively, and introduce coherence 
lengths that diverge as a power law as the fully soft state 
is approached. 

The logarithmic corrections that we predict will be 
very difficult to measure. Our universal Poisson ratios, on 
the other hand, should be observable in experiments on 
three-dimensional soft elastomers. The anomalous elas- 
ticity will be more important in 2 than in 3 dimensions. 
It would be interesting to find a two-dimensional real- 
ization of nematic elastomers, say in a crosslinked mem- 
brane confined to a substrate which inhibits height fluctu- 
ations 33] . Another two-dimensional system belonging 
the NE universality class would be a tethered nematic 
membrane that orients in a plane perpendicular to an 
external aligning field. 

Our analysis exemplifies the power and the beauty of 
the rcnormalization group. The renormalization group 
handles the rotational invariances of nematic elastomers 
in two distinct spaces quite naturally. Though the con- 
straints imposed by these invariances lead to almost baro- 
quely complicated formulae at intermediate stages of the 
analysis, our final results have a surprisingly clear and 
simple stature. 
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APPENDIX A: WARD IDENTITIES FOR 
CRITICALLY SOFT ELASTOMERS 

By virtue of the mixing invariance discussed a the end 
of Sec. lIII Al there exist several Ward identities for CSEs. 
In this Appendix we derive these identities. 

To facihtate our derivation, we introduce external 
fields via 



Tl — > Tlh = Tl — 



d'^X ^haUa + hdUd^ , (Al) 



where / d'^x — j d'^^x± J dxd- As usual, the external 
fields allow us to exploit the free energy 



F[ha,hd] = -T\nZ[ha,hd] , 
where the partition function is given by 

Z[ha,hd] = J VuaVudexpi-Tih/T) , 

as a generating function for the order parameters 

ma{-x.) = (Mq(x)) : 
mrf(x) = (wrf(x)) 



6F 
dhd{x.) 



(A2) 
(A3) 

(A4) 
(A5) 



As a consequence of the mixing invariance, the free en- 
ergy obeys the relation 

F[ha, hd] j d'^X Xahd + F[ha, h'd] , (A6) 

where h'^ — hd + 9 aha- Because the right hand side of 



d''-xxahd{x.) + d''-xmd(x.)ha{x.) 



Eq. ljA6|l is independent of the 6a, we obtain 

- 

""'Wa- 
= 0. 



(A7) 



with p 



(2) 



,(3) 



0. Another Ward identity can be 



found by taking the second functional derivative of La 
with respect to md{y) and Tnd{z). This leads in Fourier 
space to 



d 



dp^a' 



ddd 



p(i)-0 



(All) 



An identity for the four point vertex function F 



(0,4) 
dddd 



fol- 



lows in a similar manner upon taking a third order func- 
tional derivative of La with respect to, say, md{y), md(z), 
and mrf(w). The result can be stated as 



d 



' (l) dddd 
Opa 



F(",-.^)(p(i),p(2),p(3),p(4); 



— ^ add IP 



(2)^p(3)^p(4)^ 



(1,2), 



p(i)=0 

(4) „(2) „(3)\ 



odd IP ' P ' P J 



.(1,2) 
add 



(A12) 



Of course, our Ward identities should hold to arbitrary 
order in perturbation theory. At zero loop order, the 
consistency of the Ward identities can be checked without 
much effort. We carried out this check and found our 
identities confirmed. 



APPENDIX B: WARD IDENTITIES FOR 
NEMATIC ELASTOMERS I 

Here we derive the Ward identities that follow from 
the reference space rotation invariance of the NE Hamil- 
tonian H4.1|) . One of these Ward identities can be viewed 
as the origin of the soft elasticity characterizing NEs 
crosslinked in the isotropic phase. In addition, we con- 
sider semi-soft NEs where the reference space rotation 
invariance is broken by an aligning external stress. 



Now we adopt the usual strategy and take various 
derivatives. We start with 



S'La 



Srudiy) Smb{z) 



= / d^. 



S^hdi^) 



" Snidiy) Smb{z) 



dmb(z) J dmd{y)dmb[z) 

In equilibrium, the order parameters vanish for vanishing 
external fields, ha = hd = 0- Hence, we get the Ward 
identity 

- I d'x xa F^y ) (x, y , z) = Fif ) (y , z) . (A9) 
In Fourier space, identity l|A9p takes the form 



' (1) ^ ddb 
OPa 



^(A'2)(p(l),p(2),p(3) 



j<i)=0 



F 



(AlO) 



1. The soft case 

To make our arguments more intuitive we work directly 
in d = 3 dimensions. Without loss of generality we as- 
sume rotations about the x axis which we parameterize 
by a rotation angle 9. For small 9 these rotations are 
described by the orthogonal matrix 



1 
1 

-9 1 



(Bl) 



Such a rotation takes the strain u to u' — O^uQ^, where 

is the transposed of O . The relative strain w = 

— 'Ill — — 

u — u IS taken by the rotation to 

— —0 



(B2) 
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We introduce an external symmetric stress a into our 
model via 



d xaijWij , 



(B3) 



where H is the elastic energy of Eq. (|4.3II supplemented 
with the relevant bending term. Then, the free energy 



^ y = -TlnZ [g] , 
with the partition function 



(B4) 



(B5) 



is a generating function for the tensor order parameter 

SF 



my(x) = (my(x)) 



Sffij (x) 



(B6) 



Owing to the rotation invariance of the Hamiltonian (|4.3(l 
without external stress the free energy satisfies 

F\£=-9s J d'^x [(723 + ^32] +F[g:] , (B7) 

with the elements of the symmetric tensor a' given by 





= cm , 






(B8a) 


12 


= a'12 — 


0(Jl3 , 




(B8b) 


13 


= Cig + 


d<7i2 , 




(B8c) 


22 


= Cr22 - 


6 ((723 


+ (732) , 


(BSd) 


23 


= ^23 + 


((722 


- 0-33) , 


(B8e) 


33 


= Cr33 + 


((723 


+ 0-32) . 


(B8f) 



Evidently, the left hand side of the identity (|B7|I is 
independent of the rotation angle 6. Hence, the quantity 

dF r 

^ = -QQ= d''x{- sr23(x) + mi2ri3(x) 



'ni3ri2(x) - 2m23[r22(x) - r33(x) 
[m33(x) - m22(x)]r23(x)} 



(B9) 



vanishes identically. In this expression, we used the 
facts that m^- (x) and <Tij are symmetric tensors and that 
rij(x) = 2SF/Smij{x.) for i 7^ j (we use only i < j) and 
Tiilx.) — 6F/Smii{'x.) for all i (no Einstein convention). 
Now we are in the position to extract the sought after 
Ward identities simply by taking functional derivatives 
of L with respect to the order parameter. For exam- 
ple, differentiating with respect to ^23(7) and setting 
rriij = 0, we obtain the Ward identity. 



SL 



(5^23 (y) 

-s J d'*xr2323(x,y)-2[r22(y)-r33(y)] = 



(BIO) 



which in Fourier space is 

Sr2323(0 , 0)+ 2[r22(0)-r33(0)] =0. 



(Bll) 



This Ward identity is particularly important because it 
is the origin of the soft elasticity of NEs. Taking the 
appropriate derivatives of the free energy in Eq. (|4.3|l . 
it is straightforward to show that r2323(0,0) — 265, 
r33(0,0) — fli, and r22(0,0) = a2 and thus that Eq. 
(|B11|I reduces relation (|4.4a|l at 0-loop order. 

The derivation of the other Ward identities that gen- 
eralize the remaining relations stated in Eqs. 1)4. 4|l is a 
similar exercise. We restrict ourselves to stating the final 
results 



- sr2323ii(0,-p,p) - 2r22ii(-p,p) 

+ 2r33ii(-p,p) = 0, 

- sr232333(0, -p, p) - 2r2233(-P,P) 
+ 2r3333(-p, p) - r2323(-p, P) = 

- Sr231312(0,-p,p) - 2ri212(-p,p) 

2ri3i3(-p,p) = 0. 



(B12a) 
(B12b) 
(B12c) 



- Sr23231313(0,p(l',p(2)^p(3))_2r22l313(p^'\p^'\p^'^) 
+ 2r331313 (p(^\ p('\ p(') ) - 2ri22313 {p^'\ P^'\ P^'^ ) 

- ri22313(p(^\p('\p('^) = 0, (B12d) 



2. The semi-soft case 

Now we consider semi-soft NEs modeled by the elas- 
tic energy H4.1() supplemented with the explicit uniaxial 
term (|2.11|) . The aligning stress h breaks the reference 
space rotation symmetry in much the same way as the 
stress a_ utilized in Appendix lB II Thus, we basically just 
have to repeat the steps carried out in Appendix lB ll with 



Ha Ti-h.a = Hh — J d^^XUijWij 

= H-y" d^-xah.^jw,, , (B13) 



where 



C^hAj — + - {OiiOji -f Oi20j2} 0't30j3 • 

(B14) 

This analysis leads in Fourier space to the Ward identity 

sr2323(O,O)-/i = (B15) 

provided that the equations of state r22(0) — and 
r33(0) = are satisfied. At 0-loop level the iden- 
tity (|B15|) reduced to the relation between h and 65 stated 
in Eq. (j^^Uj) . 
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APPENDIX C: WARD IDENTITIES FOR 
NEMATIC ELASTOMERS II 

The mixing invariance of NEs leads to Ward identi- 
ties analogous to those for CSEs derived in Appendix IXI 
Since the derivation of the 2 sets of identities is similar 
and Appendix is fairly detailed we restrict ourselves 
here to mention differences in the derivations and to list 
results. 

Introducing external fields via the step lAll we find that 
the free energy of NEs satisfies the identity 

F[ha{^),hdi^)] ^ ~9a J d-'x Xahdix.) (CI) 
+ F[ha{Xc + OcXd, Xd), hd{Xc + OcXd, Xd) + 0aha{yi)] ■ 

Thus, the quantity La = dF/dOa obeys the equation 

La = - j dfixxahd^x) + J d'^xm.diyi)ha{x) 

f ,H , .dhJx) f . , .dhdix.) 

+ / d'^xnicix.) Xd + / d'^xrudix.) Xd 

= 0. (C2) 

By taking the appropriate functional derivatives with re- 
spect to the order parameter, we obtain the Ward iden- 
tities 

'gp{l)^ddb[P ^p(i,^o-^ab IP .P ) 

„ (2) ^ 6d IP 'P J Pa 
dPd 

„ (3) bd IP ' P ) Pa 

dpy 

j_^r^0.3)/ (1) (2) i3)^ 



(C3a) 



Jr'tL 

= ^i:,^^^(p^^^p^^O+^iy)(p(^),p(^)) 
■-42)^£"^(p^^^p'^Vi^^ 

or 



'>pd 

d 



p(i)=0 



„ (3) ^ IP ' P jPa > 

(0,4)/ (1) (2) (3) U)\ 
' (1) ^ dddd\V ; P iV : P ) 

dpi' - - 

= ^i^i^p^^^p^^^p<^v^i:..^^(p^^^p^^p 
+^il^^(p^^^p^^^p^^') 

dpy 

d 



(C3b) 



(4)^ 



(3) 

d 



dp' 



^i°i^p^^^p^^^p<'Vi^' 

\^r^dd!\p''\p''\p''')pi''- (C3c) 



APPENDIX D: ALTERNATIVE 
RENORMALIZATION SCHEME FOR NEMATIC 
ELASTOMERS 

Our renormalization scheme (|4.16(l which involves a 
re-parameterization of the coordinate Xd follows closely 
the approach developed by Grinstein and Pelcovitz |22 |. 
Though this approach has been established for more that 
2 decades now, it is not clear how the known strategies 
of proving renormalizability apply to it. In this appendix 
we briefly present an alternative renormalization scheme 
for NEs that does not entail a re-parameterization of Xd- 
We demonstrate that the alternative formulation leads 
exactly to the same results as our original approach. 

Let us revisit our Hamiltonian l|4.1()|l as a starting 
point. To reduce the number of scaling variables we set 
T = fK. This step yields 



^^-^ d'^-xj dXd{ 



•^dd 



T ~ 2f 

+ 2g VddVaa + f via +"211 ^Ib} > 



(ylud) 



(Dl) 



where w = Ci/K, g = C2/K, f = C3/K and h = C4/K. 
Note that h has no counterpart in Eq. (|4.12() and that 
we thus have an additional parameter in comparison 
to Eq. (|4.12|) . Due to this additional parameter, the 
Hamiltonian ljDl|) has an additional invariance, viz. it 
is invariant in from under the rescaling Xd P~^^'^Xd, 

Ud ^ P^/^Ud, Ua ^ P^Ua, f ^ fi^/^f, OJ ^ fJ-'iU, 

g f3~^g, f and h (3~^h (/^-rescaling). This 

scaling invariance implies that the vertex functions obey 
the scaling from 

r^^^'"'^ ({q^,g4;T,d>,.9,/,/.) = /3(*™ 

\ ^1/2^ I L II' 
"^"r^ ' 13' 13' p . 



A') 



{q±,/3 



1/2 



(D2) 



Equation ljD2|) suggest to introduce composite param- 
eters that are invariant under the /3-rescaling. Among 
various possibilities we choose 



LU = uj/h = C1/C4 , 

9 = g/h = C2IC4 , 
/ = !lh = C^/Ci . 



(D3a) 
(D3b) 
(D3c) 
(D3d) 



Note that these composite parameters are identical to the 
parameters we introduced by switching from Eq. (|4.1U|) 
to 

After this prelude we specify our alternative renormal- 
ization scheme, 



(D4a) 
(D4b) 
(D4c) 

(D4d) 

(D4e) 



T- 




- ZtT , 


UJ - 


UJ ~ 




9 - 


^9^ 


Zrj^Zgg, 


f- 




Zr^Zh / , 


h- 




Zrj. ^Zyjh . 
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With this scheme 



dxd { 



T 2T 

+ 2ZggvddVaa + Zffvl^ + 2Zhhvl^} 



Zt (viurf)' 



(D5) 



Next we derive a RGE for the vertex functions. Expfoit- 
ing as usual the independence of the bare theory of fj, and 
upon switching via t, lu, g and / to the benign parameters 
t, K, p and a we get 



[fid,. + Chdf, ^ 

X r(*^^^) ({qj 



Ptdt + + fipdp + 
,qd};t,K,p,a,h,fi) = 0, (D6) 



where 



(D7) 



and where the Wilson /3- and 7-functions are defined as 
usual. By solving the RGE (jD6|l . supplementing the 
solution with a dimensional analysis and by exploiting 
Eq. ljD2(l with the choice (3 = h for the parameter of the 
/3-rescaling we find 

r(^^'^) ({q^, 94 ; t, p, a, h, p) = h{ir'+^y^ 



1/2 



(Id 



(D8) 



Now we take a closer look at the RG flow of h. Comparing 
the renormalized Hamiltonians ljD5|) and (|4.17|) we learn 
that 



Zt Z , 



(D9) 



provided the left and the right hand side are expressed 
in terms of the same variables, here in particular t, k, p 
and a. Using Eq. ljD7p and the definition of the Wilson 
7-functions we find 



c 



-7- 



(DIO) 



From the corresponding characteristics, we consequen- 
tially obtain 



h{e) ~ z{£y 



(Dll) 



Collecting Eqs. and ||DT0|) we learn that our original 
and our alternative formulation lead to equivalent scaling 
results for the vertex functions. Of course, all results that 
follow from the scaling forms for the vertex functions, in 
particular the anomalous behavior of the elastic moduli, 
are identical for both approaches. 



APPENDIX E: CALCULATION OF FEYNMAN 
DIAGRAMS 

Here we give some details on the calculation of the 
Feynman diagrams fisted in Figs. [T] to El The first part 



of this Appendix contains 2 representative examples. All 
the remaining diagrams may then be computed by similar 
means. The second part features 2 parameter integrals. 
All 2-leg diagrams can be expressed in terms of these 
integrals. 



1. Examples 

As a first example we consider the first diagram in the 
first row of Fig.Q] For the sake of the argument we refer 
to this diagram as Add- Regarding CSEs this diagram 
stands for the mathematical formula 

^dd = qi 777^ I ka{ka + Qa) h{h + Qb) 



2 _9_ 

Id 

X Grf<i(k)G<jd(k + q) 



(El) 



where q is an external momentum running trough the di- 
agram, is an abbreviation for 1/ {27t)'^ J d'^^ k± J dkd- 
Simple power counting reveals that the superficial degree 
of divergence of Add is zero. Hence, it is sufficient to eval- 
uate the diagram at vanishing external momentum. For 
q = 0, Add reduces to 



Add ^ql\B^ M'l^ , 



(E2) 



with being a specification of the parameter integral 

m/^^ defined in Eq. ((E7|) . Exploiting the result ||ETT|) we 
obtain 



Add = qli 



51/2 



(E3) 



y4i/2 i6 7re 

as the final result for Add ■ Regarding NEs we simply have 
to replace the stemming from the vertices by {uj — gY . 

As a second example we illustrate the computation of 
the last diagram in the row line of Fig.Q] Let us call this 
diagram Cdd- In case of the CSE model it visualized the 
formula 

o2 



9 f 

Cdd = ^ j^kl [q^ • (q^ + k_L) 
X Gdd(k)Gdrf(k + q), 



(E4) 



The superficial degree of divergence of this diagram is 
2. Hence an evaluation at q = is not sufficient and 
we rather have to expand Cdd in powers of the external 

(2) 

momentum. Using the parameter integral M} ' defined 



in Eq. (|E12(I this expansion can be written as 



G, 



dd 



}(E5) 



All the other terms in this expansion turn out to be UV 
convergent. Using Eq. HE13() the final result for Cdd is 
now readily found to be 



Cdd 



^3/2 ^-£/4 



(E6) 



yl3/2 647re ■ 

Again, one simply has to replace the g^ stemming from 
the vertices by (uj — gY if one is interested in NEs. 
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2. Parameter integrals 

Many of the 2-leg diagrams can be expressed in terms 
of the parameter integral 



Jk [Bt 



[Bt + Akl + Bkl]" 



(E7) 



The calculation of this integral can be simplified by 
rescaling the momenta, 

M.W =^-V2B-(<i+'-i)/4 / ^ ^.(E8) 

Next, we employ Schwinger representation, 

(2.)^r(n) L 
X j d^^k_Lk^j^exp{-s\s.'i) J dkd exp {-skfj . (E9) 

Now the momentum integrations are straightforward. 
We obtain 



M, 



1) A-V^ g-(rf+i-l)/4p( d+f-l j 

" 2(47r)'^/2r(n)r(^) 



/ ds exp(— sBr) t 
Jo 



-1+ 



4n-l—d-l 



(ElO) 



Carrying out the remaining integration over the 
Schwinger parameter s gives finally 



M, 



X r 



(1) 



1/2 r, 2-1- 



2(47r)'^/2r(n)r(^) 

4n — I — A + S\ _ 4n-i-4+e 
T 4 



(Ell) 



In addition to Mi^^ it turns out to be convenient to 
introduce a second parameter integral, viz. 



= / — 
ik \Br 



k I 



k [Bt + Akj + Bkl]» 



(E12) 



This integral can be calculated by the same means as 
m/^^ . We obtain the result 



M, 



(2) A-3/2^-(2«-3)/2p(2i^) 

4(47r)''/2r(n)r(^) 

^ /4n — Z — 8 + e\ 4,n-l-S+e 

X r I ; 1 r 4 



(E13) 
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